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Abstract 

We compute in small temperature expansion the two-loop renormalization constants 
and the three-loop coefficient of the /9-function, that is the first non-universal term, 
for the o"- model with 0{N) invariance on the triangular lattice at N = —1. The 
partition function of the corresponding Grassmann theory is, for negative temperature, 
the generating function of unrooted forests on such a lattice, where the temperature 
acts as a chemical potential for the number of trees in the forest. To evaluate Feynman 
diagrams we extend the coordinate space method to the triangular lattice. 



1 Introduction 



Results concerning with graph theory [1-4] , that is properties of a set of points which refer 
simply to the notion of adjacency, are of interest in a variety of fields, ranging from pure 
mathematics to statistical physics and find an enormous amount of applications in natural 
sciences besides physics like in biology or in theoretical information science. 

Detailed properties of a graph can be derived from the study of the partition function of a 
q-staie Potts model [5-7] with variables defined on its sites. Indeed this function is strictly 
related with the Tutte polynomial of the graph [8-10] and, for example, the generating 
polynomial of spanning trees or unrooted forests on the graph can be recovered by taking 
the limit q ^ 0. 

A classical result in algebraic graph theory is Kirchhoff's matrix-tree theorem [11] which 
expresses the generating polynomials of spanning trees and rooted spanning forests on a given 
graph as determinants associated to the graph's Laplacian matrix. For recent applications 
see for example [12,13]. It is quite natural to rewrite these determinants as Gaussian integrals 
over Grassmann variables. 

Recently [14] it has been shown that the solution of other combinatorial problems on a 
graph can be represented in terms of Grassmann integrals, cvcnthough non-Gaussian. In 
particular, the generating polynomial of unrooted spanning forests on the graph is simply 
written adding to a Gaussian term a suitable four-fermion term. Interestingly, the same 
partition function can be obtained, order by order in perturbation theory, by considering an 
anti-ferromagnetic non-linear u-model with 0{N) invariance in the limit in which N — > — 1. 
These representations are very convenient to study the cases in which the graph is an infinite 
regular lattice, because the whole machinery of Statistical Field Theory becomes avalaible. 
For example, Renormalized Perturbation Expansion can be used, Rcnormalization Group 
notions can be applied and one sees that on two dimensional lattices these models are 
asymptotically free [15-18]. The same mapping has been used at the transition at negative 
tree fugacity which corresponds to the Potts antiferromagnetic critical point [19-21]. 

In this paper we will concentrate on the triangular lattice and, in particular, we arc interested 
in the evaluation of the so-called /3-function. We have computed the three-loop coefficient 
which is the first non-universal term, which, in contrast with the square lattice, was yet 
unknown. A direct practical relevance of this coefficient comes from a recent study of the 
zeroes in the complex plane of the partition function of the Potts model by means of the 
numerical evaluation of a transfer matrix in a strip [22]. The locum of zeroes converges to 
a pair of complex-conjugate curves with horizontal asymptote, but the convergence is very 
slow in a region of large Re(w). It turns out that the shape of this curve can be deduced 
perturbatively (in 1/w) from the expression of the /3-function, thus in the region where the 
errors are larger. 
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2 Unrooted forests 



Let G = (y, E) be a finite undirected graph with vertex set V and edge set E. Associate to 
each edge e a weight Wf,-, which can be a real or complex number or, more generally, a formal 
algebraic variable. For i 7^ j, let Wij = Wji be the sum of We over all edges e that connect i 
to j. The (weighted) Laplacian matrix L for the graph G is then defined by 



Lij 



-Wij for i=^j, 
'Ek^iWik ioTi = j. 
This is a symmetric matrix with all row and column sums equal to zero. 

Since L annihilates the vector with all entries 1, its determinant is zero. Kirchhoff 's matrix- 
tree theorem [11] and its generalizations [23-26] express determinants of square submatrices 
of L as generating polynomials of spanning trees or rooted spanning forests in G. For any set 
of vertices {ii, . . . , v} of V, let L{ii, . . . ,ir) be the matrix obtained from L by deleting the 
rows and columns ii, . . . , i^. Then Kirchhoff 's theorem states that det L{i) is independent of 
i and equals 

detL(i) = J^Il"^- (2) 

TeT eeT 

where the sum runs over all spanning trees T in G. (We recall that a subgraph of G is called 
a tree if it is connected and contains no cycles, and is called spanning if its vertex set is 
exactly V.) The i-independence of det L{i) expresses, in electrical-circuit language, that it 
is physically irrelevant which vertex i is chosen to be "ground". There are many different 
proofs of Kirchhoff 's formula ([2]); one simple proof is based on the Cauchy-Binet theorem in 
matrix theory (see e.g. [2]). 

The "principal-minors matrix-tree theorem" reads 

detL(2i, . . . , v) = ^ Yi^e , (3) 

where the sum runs over all spanning forests F in G composed of r disjoint trees, each 
of which contains exactly one of the "root" vertices ii, . . . ,ir- This theorem can easily be 
derived by applying Kirchhoff 's theorem ([2]) to the graph in which the vertices i ]^ , . . . , are 
contracted to a single vertex, while it has theorem ([2]) as a special case r = 1, through the 
bijection between unrooted spanning trees and spanning trees rooted on a given fixed vertex. 

Let us now introduce, at each vertex i E V, a pair of Grassmann variables ipi, ipi. All of these 
variables are nilpotent {ipf = tpf = 0), anticommute, and obey the usual rules for Grassmann 
integration. Writing 

'Di^,^) ■.= l[dtlj,d^,, (4) 

i£V 

we have, for any matrix A, 

J V{ij,ij) e^^^ = det A (5) 
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and more generally 



These formulae allow us to rewrite the matrix-tree theorems in Grassmann form; for instance, 
(12]) becomes 



while (131) becomes 



We 

•a=l ^ F£j^{ii,...,ir) e£F 



which is to say 



i F&T ^ i=l ' e&F 

F={Fi,...,Fi) 

This formula represents vertex-weighted spanning forests as a massive fermionic free field 
[2,27]. 

More generally, it has been shown in [14] that 

V{il),il)) exp i}Lil) + t^^ii^i + ^ WjjijjiiJiiljjiljj 



(ij) 



E (tlit\yF.\+u\EF^)\ i[we (10) 



F=(Fu...,Fe) 

where the sum runs over spanning forests F in G with components Fi, . . . ,Ff, here |VpJ and 
\Ef. \ are, respectively, the numbers of vertices and edges in the tree F^. We remark that the 
four-fermion term u U)ijipiipiipj4'j can equivalently be written, using nilpotency of the 
Grassmann variables, as —{u/2) J2i ji^ii^iLijipjipj- More interestingly, since |VfJ — \Ef,\ = 1 
for each tree F^, we can take u = —t and obtain the generating function of unrooted spanning 
forests with a weight t for each component. 



3 Relation with the lattice cr-Models. 



Recall that the A^- vector model consists of spins <Tj G M^, |crj| = 1, located at the sites 
i & V, with Boltzmann weight where 

n=-T''Y'^,,{cT,-CT,-l) (11) 
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and T is the temperature. 



Low-temperature perturbation theory is obtained by writing 



l-TTvlT'/^TTi) (12) 

with TTj G M^^^ and expanding in powers of tt. Taking into account the Jacobian, the 
Boltzmann weight is where 

n' = n + \Y.^og{i - Tn',) (13) 

i 

hj i (ij) 

When N = —1, the bosonic field tt has —2 components, and so, at least in perturbation 
theory, it can be replaced by a fermion pair ipjip if we make the substitution 

TZi ■ TZj tpi^pj - ^Jilpj . (15) 

Higher powers of tt^ vanish due to the nilpotence of the Grassmann fields, and we obtain 
the model (fTUjl if we identify 

t = -u = -T. (16) 

Note the reversed sign of the coupling: the spanning-forest model with positive weights 
(t > 0) corresponds to the antiferromagnetic A^- vector model (T < 0). 

In the case of a regular unweighted graph of order q, that is all the vertices are connected to 
other q vertices, we shall take 

{1 if z and ? are connected , 
(17) 
otherwise. 

and the corresponding Laplacian 

— 1 if i ^ j and i and j are connected 
Lij = •{ if i 7^ j and i and j are not connected (18) 
q ifi=j. 

This is the case of a regular periodic lattice in d dimensions. If we take unit lattice spacing, 
vertices connected to a given site correspond to sites at unit distance, so that, if is a lattice 
direction, / a lattice function, and x a lattice site, the lattice derivatives are defined as 

Vfc/(a;) :=f{x + ek)-f{x) (19) 
V^/(x):=/(x)-/(x-efc) (20) 
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The Laplacian can be written as 

(L/)(x):=-^V,/(x) (21) 

k=l 

and when, hke in the square and triangular lattice, q is even and to each lattice direction 
corresponds an inverse lattice direction, that is ek+q/2 — —e-ki we can restrict the sum to 
positive directions 

q/2 q/2 q/2 

-{Lf) (x) = 5] (V, - VI) fix) = Yl = E ^l^'^M (22) 

k=l k=l k=l 

q/2 

= E [/(^ + + ^(^ - ^'^) - 2/(^)] (23) 

k=l 

and, in the lattice scalar product 

(/,^) = E/(^M^) (24) 

X 

we have 

q/2 q/2 q/2 

(9: Lf) = - E(^, V^Vfe/) = E(^^^' ^^/) = E(^^^' ^If) ■ (25) 

fc=l k=l k=l 



4 The calculus of the /^-function 



We follow a procedure which already found several applications [28-32] for the square lattice. 
For a lattice theory, i.e. a theory regularized by introducing a discretization of the coordinates 
space, in principle the /9-function can be found by a direct computation on the lattice, 
which also provides a regularization. However, our lattice cr-model has a natural continuum 
counterpart, with the widely investigated action 



<S(7r,/i)=/3 j 



d^x 



(x) 



h^l - 7r2(x) 



(26) 



where we have introduced an external magnetic field h which explicitly breaks the 0{N)- 
invariance. In particular Brezin and Hikami [33] already performed the renormalization up 
to three loops in dimensional regularization. 



A general theorem of Renormalization states that the n-loop /3-function within a certain 
regularization scheme can be deduced from the knowledge of the /3-function in any other 
scheme, at the same perturbative order, and of the renormalization constants in the desired 
scheme, up to order n — 1. So, a possible procedure, which we will indeed follow in this work. 
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is to relate the /3-function on the square and triangular lattice to the continuum results of 
Brezin and Hikami via the calculation of the two renormalization constants of the non-linear 
cr-model, denoted by Z\ and Z^- 



More in detail, in our case we have to compare our lattice theory with the continuum theory 
renormalized in [33] using MS-scheme (Minimal Subtraction modified) and in dimensional 
regularization. The starting point is the relation for the n-point 1-particle- irreducible (IPI) 
correlation functions 

rLl(Pi, ■ ■ ■ ,Pn; /5, h- l/a) = Zl'^ r(;^(pi, ■ ■ ■ ZfV, Z^Z-^^l^h- y) (27) 

where a and [i are respectively the lattice spacing and the scale of renormalization for the 
continuum, while pi, . . . ,pn are the external momenta. Here we consider the lattice theory 
(denoted by subscript latt) as a regularization of the continuum theory renormalized at the 
scale l/a and we compare it with the continuum theory renormalized in the M5'-scheme 
(denoted by subscript MS) at the scale /i to determine the finite constants Zi{(3,jj,a) and 
Z2{l3,fia). Both the regularized theories satisfy a Renormalization Group equation: 

l^rg^-O; -^rLl^O; (28) 

where we added a minus sign for the lattice equation, because when a — > we are making a 
RG flux toward short distances behaviour, that has the reversed sign respect to the /i — > oo 
limit made for the continuum theory. For the lattice theory 



= -a-^r 

da 



latt 



- '^t'^iP) + (jT'-'C.fl) + 3W""{i3)) h-^ rLt , (29) 



and analogously for the MS'-theory by using W^'^'^{l3) and 7^^"^(/3) (in order to avoid con- 
fusion with the coupling costant, and in agreement with the literature on the subject, we 
denote the /3-function as W{P)). By using the condition fl27|) . we are able to join together 
the f3 and 7-function on the lattice to those in MS-scheme. Indeed we find 

W^iZ^'P) = (^Z, + W'^'\P) (30) 

7^(zr'/5) = 7"^"(/5) - (31) 

The first of them is the important relation that allows us to express the coefficients of the 
/9-function on the lattice in terms of the coefficients of the continuum theory. 

Given the /3-function for the non-linear a-model with N the number of vector components, 
we expand it in power of the coupling costant in a generic scheme of regularization 

w-'-iP) = -p - P - ^ + O(r^) ; (32) 
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the first two coefficients have not the superscript scheme because they are universal, they 
come from the calculation respectively at one and two loops (the term from order zero 
vanishes in two dimensions); explicitly they are given by 

N-2 N-2 

all the other terms are scheme-dependent; the nj^^'^'^^ coefficient is associated with l//?""''^ 
term of series expansion and correspond to a computation at {n + 1) loops. We report here 
the known results in MS'-scheme (see [33], or [29,30] for other references) 

~ 4 (27r)3 • ^^^^ 
We also expand ml/ (3 the two renormalization constants 

7(1) 7(2) 

Z.=zfH^ + ^ + 0{p-') (35) 

. X 7^^) 7(2) 
^2=4°^ + ^ + ^ + 0(r^) (36) 



With the above conventions on the series expansions, now we look at (IHUjl and we rewrite it 
as: 



^1 + i^w 



from this equation it can be seen that the coefficient of order n of the expansion of H^'"** 
(i.e. w^^^q) can be evaluated as long as one knows the coefficients of W^^ up the same order 
(i.e. , 1 ■ ■ ■ 1 w^^2) ^^"i performs the computation on the lattice of the constants Zi 
and Z2 up order — So we can argue the general result: 

W':^ = ^S::*!,oop) = F [{Wf^h={0,l,-,n-iy, {Z['^},={0,l,-,n-l}) . (38) 

For example, for the ffist scheme- dependent coefficient wi^^^, from ( 1371) we find 

wr = wo ( {Ziy - zi'A + u.iZ« + wf^ . (39) 



5 Evaluation of the constants of renormalization 



In order to obtain the perturbative expansion of the constants Zi and Z2, we use relation 
f l7r|) for the two-point function IPI. We proceed as follows: we compute r^^]^ at n — 1 loops 



^To be precise, only the constant Zi is required. The expansion for Z2 however comes out as a side result 
of the computation. 
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and, from the knowledge of at the same order, and the requirement of vahdity of flTTI) . 



we find and Z2 at n — 1 loops. 



For the continuum theory we consider the expansion 

Y%[p. /5, h- y) = -Pip' + h) + n^ip, h- /i) + 



h] /i) 



+ . . 



we report the already known two-loop results [29,30] in the case of = — 1 

h 



nSs(P' ^' =1^ (log' ^ + 8 log 7:^ - 3 + 12 (27r)2i? ) p' 
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log — + log — h 



where R is an integral defined as 



R := lim h 



dp. r dpy dq, p dq, 



h^o 2n 2n 2n J _^ 27r (p2 + h){q' + h){{p + g)2 + h) 



24 vr 



(40) 



(41a) 



(41b) 



(42) 



with '4){z) = dlogr{z)/dz, but it appears only in intermediate stages of the computation 
and cancels out in any of the results. Therefore 



Z2T^^{p,Z^'P,Z,Z,'^'h;fx) 



-1/2, 



MS 



P ip' + h) + l-{p'- h) log A + (z« - z«) p' - 



1 



^ ^og= A + 8 log A - 3 + 12 (27r)=fiV' ' 



16 TT 



log^ — + log — /l 



(43) 



1^(2) _ ^1 



(1) 



47r 



log — 

/i2 



log — 

/i 



(1) 



47r " 47r 



57r 



log — - 
^2 



47r 



57r 



6 The triangular lattice 



On a triangular lattice each site has 6 neighbours. It is convenient to introduce a redundant 
basis of three vectors e(j), as shown in figure [H such that e(j) = 0, Cj ■ ej = 1, and if i 7^ j 
then Si ■ Cj = — |. 
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Figure 1: Left: the cartesian basis and the redundant basis on the triangular lattice. Right: 
Brillouin zone in momentum space. The rhombus or the hexagon are equivalent choices, as the pairs 
of triangles denoted with b and ^ are related resp. by periodicity in ki and k2- While the hexagon 
corresponds to the direct construction of the reciprocal lattice, the rhombus is computationally 
convenient, as it is a product of one-dimensional intervals. 



Lattice sites are labelled by three integers {rii}, with x = Because of redundancy, 

a constant can be added to the n^'s without changing x, i.e. there is an equivalence relation 




(44) 



A representative of each class is chosen, for example, by fixing = 0, as 




(45) 



Remark that 




(46) 



Similarly, the conjugate quantity k 
such that Yli ki = 0. The factor | is 



= I J2i kie(i) is characterized by the three numbers k, 
introduced to have 




(47) 



As a consequence \ki\ < tt and 





so that the domain for is a hexagon of side 
of sides 27T (cfr. figure [1]). 



27r/-\/3, or equivalently a (7r/6)-angle rhombus 



Now we can introduce the Fourier transform 



f{k) for a function / on the triangular lattice 




(49) 



sites 
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which is such that 



J hexagon ^ ' 



Ih. 



(Pk 



I hexagon 

By speciahzing this general formula to the gauge (H5ll we get 



(50) 



(51) 



where we substituted = —ki — k2 and we kept into account the angle of 27r/3 between the 
vectors e(i) and 6(2) in the integration measure. Remark that the volume of the elementary 
cell generated by e(i) and 6(2) will pop out once more in the continuum limit, indeed 



sites 



X 



dki dk2 



dki 



dk2 V3 
~ ~2" 



dkrr 



-oo 27r 27r 



(52) 
(53) 



7 Tree level 



To compare with ( l26ll let us change the normalization of the Grassmann fields to get for the 
free part of the action on the triangular lattice 



-^l^tlY^ H^) [2^(a^) - + ei) - i^{x - Si)] + hti){x)i){x) \ 

sites \ i ) 



which becomes by Fourier transform 



- dk, r dk2^ 

ITT / ^M{k) 



k' + ht 



ij{k) 



where 



2sin( I 



^(2-2cosA;i) 



By using k'^ f» and fl3H|) this becomes in the continuum limit 



dk,, 



V3 J-oo 27r 2tt 
and it must be compared with the continuos expression 



(54) 



(55) 



(56) 



(57) 



(5^ 
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from which we get the identifications (see also [34]) 



A ^ ^ ht^lh (59) 



In the following it will be useful the evaluation of the integral 



in the limit of small ht- Using the relation 

COS Pi + COS P2 = 2 COS — - — cos — - — (61) 
we rewrite the denominator 

p + ht = o — Acos — - — cos — 2 cos[pi + P2) + ht (62) 



and then we make the change of variables ki — 2i±22 _ pi -pa - ^j^g Jacobian of 

the transformation is 2, but it simplifies with the factor 1/2 coming from the new area of 
integration; in fact ki and /c2 run inside the rhombus of vertices (tt, 0), (0, vr), (— tt, 0), (0, — vr), 
the Brillouin zone, which is contained twice in the square area [— tt, tt]^. So we obtain 

i(h ) = r ^ r — - f63) 

^ *^ 2n J _^ 2n 6 -4 cos ki cos k2- 2 cos{2ki) +ht ' ^ ' 

We are now able to integrate in k2 using the result 

de 1 1 



^27ra + /?cos^ - ' 



(64) 



we have 



^ht) = , ' (65) 

^ 2J(3 -cos(2A;i) + ^)2 -4cos2A;i 

dki 1 



^ 2y^(^ + 2sin2 A;i)2 - (/i^ + 9) 
2" dki 1 



27r 



2^{^-cosk,y-{ht + 9) 
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Finally, after the change cos/ci = x, we can express our integral by an elliptic integral 



(67) 



2n 



-.K 



When /i. ^ 



1 . fht 



I{h,) = -^-^ log (^^ ) + 0(/ii log /i,) (68) 



and therefore, because of (159 



We will also need the evaluation of the integral 



in the limit of small h^. Of course 



;,„,)_|_;„,)._i_ + 0(.ogA.) (71) 

and therefore 

Of course the divergent part could be obtained by going to the continuum limit 
^ dpi dp2 1 , 3 V3 p dp^ f°° dpy i i 



^ ^ — -^f^-:^ ^ ^ — : — (73) 

/ 27r 2n [p^ + h^f 2 2 1^ 27r J.^ 27r f [p2 + 4^- ^ ^ 



1 TT 



^From 3.148.2 of [35] 



dx = =F{p, r) (66) 

d ^y {a — x){b — x){c ~ x){x — d) y'{a — c){b — d) 



with a > b > c> u > d and /3 = arcsin y ^ ~ y (a-c)(b-d) • -'--'^ ~ ^'''2^^ + V^t + 9, 5 = 

F{P,r) = rf* , '^^ is the elliptic integral of the second kind, and if /3 = , F{^,r) — K{r) is called 

^ -y/ 1 — sin^ ^? 

the complete integral. 
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> X > 




Figure 2: The Feynman diagrams for the two-point function at order 1. 



Analogously 



lim ht 

ht^O 



dpi dp2 dqi dq2 1 
- 2vr 271 2iT 271 + ht){q^ + ht){f+^q' + ht) 



~ lim h , 



dp^ dpy dq^ dqy 



1 



R 



3 J 2ti 2ti 27r 2ti (p2 + h){q^ + h){{p + qf + /i) 

where i? was defined in fH2l). 



(74) 



8 One-loop diagrams 



The interaction terms on the triangular lattice are 



- ^ / V-lQ + fc)^(g) Hp - kmP) (75) 

^ Jp,q,k 



where we introduce the shorthand 

dfci r dk 



2. 2. ■ 

We wish to compute the IPI two-point function. At one loop, two graphs contribute (fig. [2]). 
On the triangular lattice, by defining 

A(A;) := + ht (77) 

we get 



i-fi-i + hti+ If (1 - hti) 

o 
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By going to the continuum limit, in the hmit of small magnetic field 

2 r 3 p2 3 1 ha^} 

+-^[p^-/.]log^ (79c) 



2^3 47r J ^ 48 

By comparing the two expressions we obtain the one-loop result 



which, of course, result to be independent from the magnetic field. 



9 Two-loop diagrams 



The diagrams at second order are the four ones shown in figure [3l As we expected these 
are the same Feynman diagrams that appear at the second order of perturbative expansion 
of the (T-model [29]. According to the Feynman rules we have to add a minus sign to the 
diagrams A and C: for the first one since it has a mass insertion, for the second since it has 
a loop. So that the expression of the second order contribution of self-energy is 

ni = A- B + C- D (82) 

with 

^ = A,, /\{p + q)A{k)A{k + q) ^^^^ 

n = /" P-^g /ofjN 
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Figure 3: The Feynman diagrams for the two-point function at second order. On the top left 
corners, we report the identificative letters. 



The first two diagrams are easy to evaluate exactly in terms of I and l2- We find 



ll + I^ + hth 



+ / - hth 



B =f 



-J + I2 + Ihth - 2htl2l + 4 

6b 



2 3 
+ 21 -\- 3htP - hth + -hj ^ 



+ 



^hlh + 2h1hl 



6 " ^ 2 

The diagrams C and D are more involved. First of all remark that 

[/\{p + k + q) - /\{k)] 



D-C 



k'^[p + k + q -k'^] 
,^^A{q)A{k + q)A{p + k) " j,^^ 

A{k) [A{p + k + q) - A{k)] 
fc,, Aiq)Aik + q)Aip + k) 

A{k) [A{p + r) - A{k)] 
kr A{r - k)A{r)A{p + k) 



A{q)A{k + q)A{p + k) 



(87) 
(88) 
(89) 

(90a) 
(90b) 
(90c) 



where in ( I90bp we neglect terms of higher order in the small-/i expansion. We are interested 
in the first terms of the Taylor expansion for small external momentum. We get 



D 



A(r - A;) 



sin Pi sin ki sin pj 



sin ki 



sm Ti 



A{k) A(r) 



A(r) A(r) 



p 



1 

2 AM 



A{k)A{r) 



sin Pi sin k, 



(91) 
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We easily get 



1 



Ik,r A(r - k) 
_ r A{k) 
Jk,r A(r - /c)A(r) 

r 9 

A,, A(r - A;)A(r) 

■ 



= I 



6 3 



— P I 



2 X,, A(r - k)A{r) 

We have still to compute (changing r into — r) 

. . . /" sin fci 

4 > sm »j sm pj / — twttt 

7,,, A(r + k)A{k) 

above is 

sin Pj sin 

i.i 



The tensor form of the expression 



sin kj sin r j sin kj 
A{k) ^ 'Airj ~ A(r) 



(92) 
(93) 
(94) 
(95) 

(96) 
(97) 



with Ajj symmetric under the exchange of i with j, and permutation of indices 1,2, 
in general Aij = a + b6ij, which substituted into the previous expression gives a(X^j 
b^- sin^ Pi ~ bp^ + 0{p'^) because we have that Y^^Pi = 0. Therefore we need 
coefficient b which can be computed, for example, as 



sin 



Then we get 



f sin ki — sin k^ 



f (sin ki — sin k^) sin ki f (sin /ci — sin k^ 

^ A(r + k)A^{k) ^ ^ A^ 

= I [1(7 - /i,/2) - ^ 



oiii /ci sin ri sin ki 
A{k) ^ A(r) ~ A(r) 

(sin /ci — sin /cs) sin /ci 



(98) 



and 



f sin^ /ci 

A,, A(r + A;) A(A;) A(r) 

/■ sin /ci sin 
A,, A(r + ^) A(A;) A(r) 

/■ sin Aci sin ri 
^ A(r + A;)A(fc)A(r) 

/" sin sin ri 
A,, A(r + *;)A(A;)A(r) 



3 V6 2V37r, 

~"6 ^ Ws^^Ts ■ 



12 ^ sVStt. 

_ ^ _ G 
~ 9" ~ T 
G _ J_ 
8" ~ 144 



6 ^ Vl2 4y37r. ^ 



R G L 
18 ^ ¥ ^ 24 



P I R G K L 1 
12 ~ 8^3 TT ~36^16~16~48^288 



(99) 
(100) 
(101) 
(102) 
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with 



G 



K 



ki + ri [A(r + k) - A{k) - A(r)] 
k,r A2(r + A;)A(A;)A(r) 

ki k2k'ifiki + ri 
,,,A(r + A;)A(fc)A(r) 

ki + ri ki Ti 
A(r + fc)A(A;)A(r) 



(103) 
(104) 
(105) 



So finally we found: 
and in conclusion 

ni=f 



2V3 7r 



+ + 



3G + ir + L 



1 ■ 
72 



1 
6 



/ + 9 + /it/' - ^/ 



3 



(106) 



21 1 

+ + 



VStt 72 24v^7r 
1 



-2r 



+ R + 
1 



3G' + is: + L 



(107) 



6 2^3 tt/ 24^3 tt. 
By comparing the two expressions ( l43l) and (11071) we obtain the two-loop result 



(2) 



9 



2 /^'a' 



^(2) 



167r2"'° 48 ' 47r 

5 , n /i^a2 1 
log' ^ + 



+ 



87r' 



log 



/i2a2 3(3G + ii' + L) 



48 



+ 



+ 



167r2 8 



+ 1 (108) 



167r2^^^ 48 ' 4v^7r 
and the three-loop result by 



log 



Wo 



latt 



levr 87r2 



(109) 



IGtt^ 



9{3G + K + L) 
8^ ' 



(110) 



By application of the coordinate-space method by Liischer and Weisz [36] suitably modified 
for the triangular lattice (see appendix we have obtained the numerical determinations 



G = -0.025786368 
K = -0.007632210 
L = 0.035410394 

with errors smaller than the quoted digits, from which we recover the value 



(111a) 
;illb) 
(111c) 



w 



latt 



-0.01375000819, 



(112) 
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10 A direct application 



The determination of the coefficient w'f^* can be used, as shown in Ref. [22], to recover, for 
example, the phase boundary in the plane of complex temperature for the g-state Potts model 
in the limit g — > 0. This separatrix is, indeed, a special renormalization-group flow curve. 
If we call X and ?/, respectively, the real and imaginary part of the complex temperature we 
must have therefore that 

y{x) = y, + ^ + ^ + (113) 



X x^ 



where 



and I/O "was numerically estimated to be 

I/O ^ 0.394 ± 0.004 (115) 
For numerical purposes in [22] a variant parametrization is followed, that is 



y{x) = yo exp 



Bi B2 B; 



H = \ h 



(116) 



X — (^x — aoj {X — 
where comparison with (11131) in the limit of large x gives the relations 

B,=A, = --^ B2=A2-^-aoA,. (117) 

The parameter ao, and Ai and therefore Bi with i > 2 were not known. In [22] the authors 
decided to truncate (11161) by setting _Bj = for i > 3 and try estimated and B2 by this 
ansatz by imposing the value of the function and its derivative on the last known numerical 
point, that is ?/(0.0198) = 0.23 and ?/'(0.0198) = 0.369003. They estimated 

ao ~ -0.550842 (118) 
B2 ^ -0.122843. (119) 

From our calculation we get an evaluation of B2 

^0.00959932 (120) 

B2 ^ 0.0053776 + "V (121) 
2nV3 

so that we can use the strategy just discussed to derive B^ in addition to ao and B2. We 
obtain 

ao ^ - 0.778527 (122) 

B2^ - 0.066160 (123) 

53 ^ - 0.162495. (124) 

The curve resulting from this numerical values does not differ substantially from the old one 
as can be seen in Fig. HI This gives more confidence on the method and results in [22]. 
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0.4 



5 



Im(l/t) 




0.003 




0.2 



0.002 



0.001 











12 3 4 



-0.2 



-0.4 



Re(l/t) 



-0.2 







0.2 



0.4 



Figure 4: Phase boundaries for infinite strips of the triangular lattice. Numerical values from 
different lattice widths L, from 2 to 9, on gray-tone curves from left to right. Black dots reproduce 
the extrapolated L — > co limiting curve in the region of negative Re The black dotted-dashed 

curve and the continuous black curve (almost indistinguishable), in the region of positive Re (1/t), 
are respectively the old and new curves from the ansatz of equation ()116p . In the magnification on 
the right, we plot the discrepancy 6 between the two curves along the Im(l/t) axis, as a function 
of Re as it should, it vanishes with its first derivative, at the numerical point 0.0198 used for 

the extrapolation, and vanishes asymptotically because the same estimate of the asymptote yo is 
used; all in between, it remains of order 10"'^. 
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A Liischer-Weisz method for evaluation of lattice in- 
tegrals 

In the evaluation of two-dimensional lattice integrals, we used the coordinate method il- 
lustrated in the paper by Liischer and Weisz [36], and specialized to two dimensions by 
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Dong-Shin Shin [31,37], although also the momenta method proposed in the appendix C 
of [38] could have been used. 

The main idea is the use of some basic relations for the free propagator in coordinate space 
(the defining Laplacian equation and a set of relations due to Vohwinkel) , in order to find a 
recursion which, starting from the values in a certain number of sites neighbouring the origin 
(the fundamental lattice integrals) , allows to find the whole set of free propagators in lattice 
sites in a large radius R, in a time which scales polynomially with R. As a side result, it 
gives a simple proof of the fact that all these values are linear combinations with rational 
coefficients of the fundamental lattice integrals. 

Generalization of the procedure to the triangular lattice is not straightforward, and involves 
some delicate points. Some of them are: 

• In the redundant set of variables (pi,P2,P3), the constraint J2iPi — does not allow 

for derivatives in a single variable: one should either perform linear combinations of 
derivatives where the sum of coefficients is zero (for example, (Vi — ^3)fiPi,P2,P3)), 
or equivalently, perform derivation within a non-redundant choice of variables (for 
example, Vi/(pi,p2, -pi -P2))- 

• because of this fact, the Vohwinkel relations involve a larger number of terms, and thus 
it is more difficult to manipulate them in order to have a recursion relation. It will 
turn out that a larger strip is required for the first x-axis recursion. 

• For Ge{x) at values of £ larger than 1, the choice of subtraction is now not anymore 
easily deduced by the Taylor expansion of the exponential and the requirement of 
periodicity. Now we also have the requirement of gauge-invariance under x ^ x + 
m(l,l,l), which forces the application of "hat" factors only to combinations of pi 
where the sum of coefficients is zero. 

Coming back to the point, the free subtracted propagator 




ik' 



(125) 



statisfies the Laplace equation 



AG{x) 



5(2) (a;) . 



(126) 



with the lattice operators 



3 



3 




(127) 



i=l 



i=l 



and we have 



^(0,0) = 



G(1,0) 



1 



(128) 
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For the triangular function H{x) defined as 

H{x) = I e'^'^lYiP (129) 



a set of Vohwinkel relations holds 

G{x + /i) - G{x -T) - G{x + v)+ G{x -v) = (x^ - Xy)H{x) (130) 
but only two of them (e.g. (/x, v) = (1, 2) or (1, 3)) are independent. 

Using the previous equations and the Laplace equation fll26p . we are able to eliminate H{x), 
and write a recursion relation. The one we find on a width-2 strip along the x axis is given 
by the set of equations 

= -6G'(x, 1) + ^^((0;, 1) ±/i) (131) 

= -6 G{x, 0) + 2{G{x +1,1) + G{x, 1)) + G{x +1,0) + G{x - 1, 0) (132) 
= {G{x - 1, 1) - G{x + 1, 1)) + x{G{x, 2) - G{x, 0)) 

+ (x-l)(G'(x + l,2)-G(x-l,0)) ^ ' 

which must be solved with respect to G{x+1, a), with a = 0, 1, 2, in order to have a consistent 
recursion. A new fundamental integral is required. A choice could be G{2, 1), which is valued 

G(2,l) = l-^ (134) 

In a similar fashion, given the values of G{x) on the width-2 strip, the function can be 
determined in the whole plane (a sector with x = {n,m), with n > 2m > is sufficient, 
because of symmetry). The Laplacian equation alone is enough to fulfill this task. So we 
conclude that at all values of x the function G{x) is in the set Q + ^Q. 



The integrals of the form 

Kf"-^ = / ^ (135) 
J p P 

which involve G{x) only on the real axis, are easily computed, the first values being 
Kt = ~U'-^ Kf^m-?^ , _ _448 ^ 288V3^ 

3 TT TT 3 TT 

The next ingredient we need in order to calculate all the triangular-lattice quantity arising 
from our diagrammatics is the two-propagator function in coordinate space. It turns out 
that the proper subtraction is the following 



kx 



1 + if (^^ - 2A;|)(xi - X2)^ + cyclics 



G.(x) = / ^ (137) 
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The triangular-lattice Laplacian relation still reads 

- AGsCt) = G{x) (138) 
while the Vohwinkel relation, still for (/i, z/) = (1, 2) or (1, 3), is 

G2(x + Ji) - G2(x - Ji) - G2(x + u) + G2(x - 9) = -(xn - Xi,) (g(x) H ]=—] (139) 

V 4v3 7r/ 

(remark the presence of the corrective contribution l/(4-\/3vr) due to regularization). At the 
aim of building the recursion, also in this case it turns out that, as the "support" of the 
relations is identical to the one of the triagular-lattice G{x) case, the independent lattice 
integrals still must be the ones located at the points x G {(0, 0), (1, 0), (2, 1)}. The first two 
vanish because of the subtraction, while the last one is computed analytically, with the result 

G2(2,l) = -^ (140) 

and thus, as it is again a rational times ^/S/tt, still the function G2{x) at a generic point is 
in the set Q + ^Q. 

The two-propagator analogues of the quantities Kf"^ are the integrals of the form 

K2n) _ / Pi 



They still involve G2{x) only on the real axis, and thus are easily computed, the first values 
being 

= /f| = -4 + ?^ = (142) 

3 37r TT 3 77 

We need also the lattice sums 

G=Y^[{Vl-V,fG{x)][G{x)f (143) 

K = 2J2 [(VI - Vi) (V; + V^) G{x)] [{V; + V3) G{x)] G{x) 



+ E [(Vt - (V; + V2) (V* + V3) G{x)] [G{x)Y 



(144) 



L=Y. [(Vt - Vi) G{x)f (145) 

xez^ 

At this aim we need the full strength of coordinate method: we evaluate the subtracted 
propagators G{x) and G2{x) on lattice points up to a given hexagon of side r (~ 25), exactly 
in terms of rationals, in negligible computational time {0{r^)), from which we deduce the 
(9(r^) largest terms in the sums above, while the remaining contribution is estimated from the 
large-distance behaviour of the integrands. The numerical results are reported in equations 
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B Details on two loop lattice integrals 



By using the identity 

2 [sin a + sin + sin 7] = —a /3 7, (146) 

valid when a + /3 + 7 = 0, in the cases (a, /9, 7) = (fci, /c2, /C3) or (fcj, r^, —ki — Vi), we easily 
get 



/■ fei k2 h ^ _2 r A _ 

A,,A(r + A;)A(^)A(r) ^ J 



J,^,A{r + k)A{k)A{r) J,^, A{r + k)A{k)A{r) \ tt J 6 

f EijSinfejSinr,- /■ hnk2r2hf3 f 2 V3\ 1 

A,,A(r + fc)A(fc)A(r) A,, A(r + fc) A(fc) A(r) 1^ tt 7+'''^ 2 + 



[sin /ci + sinri — sin (/ci + ri)]^ /" /ci + ri ki f^^ 



12 



A(r + fc)A(fc)A(r) A^,, A(r + k)A{k)A{r) 

nj=i 2 [sin + sin - sin {ki + n)] /" kifih[^ik2f2h^^2 L 3K 



,k,r A{r + k)A{k)A{r) Jk,r A{r + k)A{k)A{r) 2 2 

We also see that 

-2 ^ [sin ki + sin ^ - sin - {ki + n)] = - ^ fcjfjA;, + n 

= /ciA;2/i;3 + fif^2r3 - ki + rik2 + r2k^ + 

therefore 

/" hfiki + ri E- fc^f^/cj + ri _ r -kik2h Y.i hrjkj + Tj _ 
A,, A(r + A;)A(A;)A(r) "A,, A(r + A;)A(A;)A(r) " ^ ^ 

computed either using the first two lines or the last two lines of the previous block of 
identities. 
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